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Abstract — This paper presents the use of Inverse Free 
Krylov Subspace Algorithm (IFKSA) with Estimation of Signal 
Parameter via Rotational Invariance Technique (ESPRIT), for 
Direction -of -Arrival (DOA) estimation using a Sparse Linear 
Array (SLA) of antenna elements. The use of SLA reduces the 
hardware requirements and the production cost. The SLA 
generates the sparse data, the pure signal subspace is first 
decomposed from the received corrupted signal space using 
IFKSA and later, ESPRIT is used to estimate the DOAs. 
Computer simulations are used to evaluate the performance 
of the proposed algorithm. 

Index Terms — Krylov subspace, ESPRIT, DOA, Sparse linear 
array 

I. Introduction 

Estimating the Direction-of- Arrival (DOA) of the far-field 
narrowband signals has been of research interest in array 
signal processing literature [1]. The applications of the DOA 
estimation is found in RADAR, SONAR and other Wireless 
communication systems in the areas like localization, tracking, 
navigation and surveillance. There are number of DOA 
estimation algorithms like Multiple Signal Classification 
(MUSIC) [2], Estimation of Signal Parameters via Rotational 
Invariance Technique (ESPRIT) [3] and their derivatives 
[4] [5] [6]. These algorithms rely on the uniform structure of 
the data, which is obtained from the underlying Uniform 
Linear Array (ULA) of antennas (sensors). ULA consists of a 
number of elements with inter- element spacing equal to half 
the wavelength of the impinging signal source. The size of 
the ULA increases the hardware requirement. The same 
performance with slight modification in the Mean Square 
Error (MSE) can also be obtained by using the Sparse Linear 
Array (SLA) of antennas or the Nonuniform Linear Array 
(NULA) of antennas[7][8]. In the SLAs, the inter-element 
spacing need not be maintained and the reduced number of 
elements further reduces the hardware requirement and the 
production cost. Moreover, the SLA gives the same 
performance as the ULA, with reduced number of elements. 
Furthermore, failure of a few elements at random locations 
will also result in SLA. 

Among the various high resolution methods for DOA 
estimation, subspace based methods are most popular and 
powerful. The popularity is due to its strong mathematical 
model to illustrate the underlying data model and it can with 
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stand the perturbations in the data. Subspace methods for 
DOA estimation searches for the steering vector associated 
with the directions of the signals of interest that are orthogonal 
to the noise subspace and are contained in the signal 
subspace. Once the signal subspace is extracted the DOAs 
are estimated. The decomposition is performed using the 
Eigen Value Decomposition (EVD) of the estimated received 
signal correlation matrix. The MUSIC and ESPRIT are the 
popular subspace based DOA estimation algorithms. MUSIC 
algorithm is a spectral search algorithm and requires the 
knowledge of the array manifold for stringent array calibration 
requirement. This is normally an expensive and time 
consuming task. Furthermore, the spectral based methods 
require exhaustive search through the steering vector to find 
the location of the power spectral peaks and estimate the 
DOAs. ESPRIT overcomes these problems by exploiting the 
shift invariance property of the array. The algorithm reduces 
computational and the storage requirement. Unlike MUSIC, 
ESPRIT does not require the knowledge of the array manifold 
for stringent array calibration. There are number of variants 
and modification of ESPRIT algorithm. ESPRIT algorithm is 
also extended for sparse linear antenna arrays or non-linear 
antenna arrays. For nonlinear arrays the aperture extension 
and disambiguation is achieved by configuring the array 
geometry as dual size spatial invariance array geometry [9] or 
by representing the array as Virtual ULA, and using the 
Expectation-Maximization algorithm [10]. However, the 
subspace algorithms are heavily dependent on the structure 
of the correlation matrix and are unsuitable to handle sensor 
failures. 

In this paper, we extend the conventional ESPRIT to 
estimate the DOAs using the SLA. The basic idea of any 
subspace based approach is; decomposing the received 
corrupted signal space in to pure signal subspace and noise 
subspace. This is generally achieved using the Eigenvalue 
decomposition. However, the sparse data output from the 
SLA cannot be directly decomposed, because it deteriorates 
the performance of the subspace based approaches. We 
here, use the Krylov subspace based methods [11] to 
decompose the received corrupted sparse signal space and 
later use ESPRIT to estimate the DOAs. The resulting 
algorithm is called as IFKSA-ESPRIT DOA estimation 
algorithm. There are number of Krylov subspace based 
techniques. The iterative methods like Lanczos algorithm, 
Arnoldi algorithm and Jacobi-Davidson algorithm are very 
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popular among them [12]. The methods are very sensitive to 
perturbations and are shift-and-invert procedures leading to 
increased computational complexity. Inverse Free Krylov 
Subspace algorithm proposed by Golub and Ye [13], is an 
attractive eigenvalue computation algorithm. IFKS A iteratively 
improves the approximate eigen pair, using either Lanczos or 
the Arnoldi iterations at each step through Rayleigh-Ritz 
projection procedure [11]. The algorithm is a very attractive 
due to the following reasons; first, the technique can be used 
to find any number of smallest eigenvalues (Largest can also 
be calculated), and second, the algorithm is less sensitive to 
perturbations. The performance of the algorithm is evaluated 
for various sparse antenna array geometries. 
The rest of the paper is organized as follows. In the following 
section the signal model is discussed. In section III, the 
proposed DOA estimation algorithm is presented. Section IV, 
discusses the simulation results and finally, conclusions are 
drawn in section V. 

II. Signal Model 

In this section signal model is discussed, beginning with the 
ULA and later we discuss for the SLA. The DOA estimation 
problem is to estimate the directions of plane wave incident 
on the antenna array. The problem can be looked as parameter 
estimation. We here mainly introduce the model of a DOA 
estimator. Consider an M - element uniformly spaced linear 
array. The array elements are equally spaced by a distance J, 
and a plane wave arrives at the array from a direction & off 
the array broadside. The angle & is called the direction-of- 
arrival (DOA) or angle-of-arrival (AO A) of the received signal, 
and is measured clockwise from the broadside of the array. 

Let N narrowband signals s L (£}, s 2 (0, -• >s N (t) all 
centered around a known frequency, impinging on the array 
with a DOA ft . I = 1,2, ■ ■ ■ N . The received signal at the 
array is a superposition of all the impinging signal and noise. 
Therefore, the input data vector may be expressed as 



Where , 

a (0.}, is the steering vector of the array in the direction 
& t . To further simplify the notation, we write (2) as 

X = A(Bjs-h W (3) 

Where , X is the array output matrix of size if X K , A is 
the complete steering matrix of size M x JV function of the 
DOA vector B , is signal vector of size and is the noise 
vector of size M X K . Here, K is the number of snapshots. 
Here, (3) represents the most commonly used narrowband 
input data model. 

We are primarily interested in a SLA, where, some elements at 
some locations are missing. An example of a SLA is shown in 
Fig. 1, where elements at locations 10, 7 and 3 are missing, 
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Figure. 1. Example of a SLA. Elements at Location 3, 7 and 10 are 
Missing. 

and the remaining elements are feeding the output.The % 
output of a SLA will be represented as a sparse data given 
as 

X = A(8)sfW (4) 



Where, X^ x ^ is the sparse data generated by SLA of 
Af elements, B) ^ yJl r * s tne steering matrix of the SLA, 
and Wft xx is the noise. 

The estimated correlation matrix of (4) of the sparse 
represented data is written as 



h! = l 



(5) 



Therefore, the problem is to decompose the correlation matrix 
in (5) to obtain the signal subspace, which is further 
processed using the ESPRIT, to estimate the DOAs. 

III. Ifksa-Esprit Algorithm 

The IFKSA-ESPRIT algorithm for estimating the DOAs from 
a faulty ULA is given in Table I. It consists of two steps. Step 
1, is decomposing the correlation matrix given in (5) using 
the IFKSA technique. In the second step, the conventional 
ESPRIT algorithm is applied to estimate the DOAs.The 
IFKSA algorithm for finding the approximate N largest 

eigenvalues and the eigenvectors Qf.v)°f me correlation 
matrix^ given in (5) is described as follows. The IFKSA 
algorithm computes the smallest eigenvalues of k • To 
compute the largest eigenvalues replace by r and - R 
compute its smallest eigenvalues and reverse the sign to 
obtain the largest eigenvalues of [13]. 

Given an initial approximation, the goal is to improve the 
initial approximation through the Rayleigh-Ritz orthogonal 
projection on a certain subspace, i.e., by minimizing the 
Rayleigh quotient on that subspace. 



(6) 



The IFKSA is a technique for computing the new approximate 
eigenvectors vv +1 .The technique starts with an approximate 
eigenvectors v -. and constructs a new 
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Table I. Inverse Free Preconditioned Krylov Subspace - Esprit 
Algorithm (IFKS-ESPRIT) 

Input R: 

Output: The estimated DOAs 

Step 1: Compute the Eigen decomposition of the array 
covariance matrix using the Inverse Free preconditioned 
Krylov subspace (IFKS) algorithm 
IFKS Algorithm ( ): 

Step a: Set n£ 1., Rand an initial approximation wifh 

Nf=i 

Step b : For it = 0,1, ■ ■ -, unti Z l l onverg enze 

Construct a basis = [u^ti^ ■■-..ti.J for the 
subspace 

JT =^n(^,(a- P]E )^,(R 

- VJ 

Step c: Foim E m = U m 7 (R - pju m 
Step d: ComputE th» CftjWj^Knpaiis of K^. and 
select the desired ones. 

End; 

Step 2 : E stimate the DOA u sins E SPRIT 
ESPRIT (): 

Step a: Form the siaial subspace V, usuie ftie 

=iE=-nvei:tors obtained from step 1 
Step b : Form the matrices V lfi and V, 
Step c: Estimafe feDOA bdRthe [V v V.J . 

End; 

□ 

approximate by Rayleigh-Ritz projection of on to the 
Krylov subspace. 

JT 1 = span (vfe, (R - pjifc, (r - p J ^ ■ -- , (r - pOV) (7) 

where ^ is the fixed parameter, generally taken between 
1 < n < 100 [15]. 

The projection is carried out by constructing a matrix of 
basis vector for 3€ then forming and solving the 
projection problem for . The process is repeated for the 
required number of eigenpairs. This iteration forms the outer 
iterations. In the outer iterations, using Q n : 

* ra =U m T (*- Pk )U m (S) 

is formed. And the eigenpairs Ql^ w t ~) for R is computed. 
Then the new approximate of the eigenvector is 

= U^il (9) 

and correspondingly, the Rayleigh-Ritz quotient 

pK + i=Pft+A (10) 
is a new approximate eigenvalue. Now, to construct the or- 
thonormal basis vector y ^ he Lanczos algorithm [1 1] is used 
as inner iteration, the algorithm is given in Table 2. The 
Lanczos algorithm starts with an initial vector v D One matrix 
vector multiplication g= C-^u^ performed in every iteration, 
where is the residual. To make sure that + 1 is orthogonal 
to u j -iand U j and extra orthogonalisation against these two 
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vectors are done in further steps shown in Table II. 

The steps of the IFKSA are repeated until convergence. 
Once the signal subspace is computed from the N 
eigenvectors v k 5 the DOAs are estimated using 

sparse Datn 
X 

^> 



Compute ifte ConeJatiofi 
Malrix of X 
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IFKSA [o deawnpose- ihe 
spars* data and obiahn the 
pure signal subspace 
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□OAS 
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Figure2: Implementation of the IFKSA-ESPRIT Algorithm. 
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Figure 3: Flowchart of the Proposed Algorithm 

conventional ESPRIT method. The implementation of the 
proposed algorithm is illustrated in Fig. 2. The flowchart of 
the proposed IFKSA-ESPRIT algorithm is shown in 
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Fig. 3. The flowchart starts with initializing the maximum 
number of iterations, for every iteration the estimated 
correlation matrix of the sparse represented data for K number 
of snapshots is computed. From the computed estimated 
correlation matrix, the large N eigenvectors are calculated 
from the IFKS A algorithm. Using these N 

Table II. Lanczos Process For Constructing The Orthonormal Basis 



Input C t = (r — pjj, an approximate siEaLvaatDf: \\ 
•Qr^i3l{U = [ii y ii 3J - uj 

for l = 0: (n - l) 

S = C t U: 

«j = ujg 

II: t 

^ = /F; 
End: 



Table. III. Summary Of The Ifksa-Esprit Algorithm With Computational 
Complexity 



Algorithm 


Computational 
Complexity 


R = i[xx«] 




Calculate [ *iE*nvalu*s . ?iE*nv actors] u sinE IFKS A 




method 


+ 1} 


SalKt X larEast ?i|?n pair to obtain the siEnal 


subspaceV, 




Fonn V rf by dalatinE ihe first column of V a 




Form V al by delating tha last column ofV, 




Compute EiEam-aluas 




= eiE,(V rf jiinv(V al )) 




Calculate DOAs 





larger eigenvectors the conventional ESPRIT is executed to 
estimate the required DOAs. The results are computed and 
stored for each iteration for evaluation, until the maximum 
number of iterations. The summary of the algorithm with the 
computational complexity is presented in Table III. 



IV Simulation Results 

In this section, we examine the performance of the 
proposed IFKSA-ESPRIT algorithm. The signal model 
consists of two equal magnitude complex analytical signal 

impinging on a 20 element ULA from direction [10°, 4°. 

The interelement spacing in ULA is equal to half the 
wavelength of the signal. The received signal is noncoherent 
and corrupted due to Additive White Gaussian Noise 
(AWGN). All examples assume 200 snapshots and 100 
simulation trails are conducted. The simulations are performed 
using Matlab 7, on a Intel core 2 Duo processor installed 
with Windows Xp SP2 operating system. 

The signal is assumed as a complex exponential sequences, 
given by 

s(t) = exp (jtpiOl t = 12, -"K 
Where is the uniform random number distributed 
between - n: to n: The Signal-to-Noise Ratio (SNR) is defined 
as 
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SNR =lOto0 1D -^ 

where, tr£ and are signal power and noise power 
respectively. In the first example, the performance of the 
proposed IFKSA-ESPRIT algorithm is compared with the 
ESPRIT algorithm and Cramer Rao Lower Bound (CRLB) [14] . 
The Root Mean Square Error (RMSE) of the estimated DOAs 
are calculated and plotted with respect to various Signal-to- 
Noise Ratio (SNR). The RMSE is calculated using the 
equation, 

RMSE =tJe{(& -§Y} 

From Fig. 4, it can be seen that for the case of ULA the 
IFKSA-ESPRIT is showing improved performance when 
compared with the ESPRIT algorithm and is close to CRLB at 
low SNR. In the next example, the performance of the 
proposed algorithm is evaluated for SLA. Various SLA 
geometries assumed for the simulation is listed in Table. IV. 
The numbers in the table gives the information on the position 
of the elements feeding the output and the missing numbers 
means the corresponding element is missing. The estimated 
DOAs are also plotted on histogram with respect to number 
of trails. The histogram plots give the information on the 
accuracy of the estimator to estimate the DOAs. From the 
RMSE vs. SNR plot for all the four cases is shown in Fig. 5, it 
can be observed that except for case 2, in all other cases the 
algorithm is able to estimate the DOAs. It can be observed 
from Table. IV, for second case the number of elements are 
very less, resulting in performance degradation. The 
histogram for all the three cases (See Fig. 6 - Fig. 9) show that 
for the second case of SLA the accuracy is poor, see Fig. 7. 

Table IV. Elements Considered To Form Sla In Simulation 



ELEMENTS CONSIDERED TO FORM SLA IN SIMULATION 



SLA 


Positions 


SLaI 


[1.2.14.:-. 6.15. lb". r.lB. 15.50] 


SLA2 


[1.2,3.4.544.17.20] 


SLA3 


[I.2.3.4.5.E. 12.1I5.1B.20] 


SLA4 


[1^^7,10,12,13.14,1 6.1 7.19.20] 


10°, ■ ■ ■ , 




10 1ft 20 

Figure 4: RMSE vs. SNR plot of IFKSA-ESPRIT, ESPRIT and 
CRLB 



ACEEE 



Figure 9:Histogram of the estimated DOAs fro case 4. Assumed DOAs of two signals are [1CT.4-. SNR = 10 dB 

ACEEE Int. J. on Communications, Vol. 03, No. 02, Nov 2012 




II 5 10 IS 20 55 X 3& *3 

Figure. 5. RMSE vs. SNR plot of estimated DOAs for all the 
four cases of SLA. Assumed DOAs of two signals are[lO* 45*] 
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Figure. 6. Histogram of the estimated DOAs fro case 1. 
Assumed DOAs of two signals are [icp, 45 3 ] SNR = 10 dB 
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Figure. 9. Histogram of the estimated DOAs fro case 4. Assumed 
DOAs of two signals are [10* 45 3 ] SNR = 10 dB 

V. Conclusion 

In this paper, we proposed IFKS A-ESPRIT DOA estimation 
algorithm which uses the Sparse Liner Array of antenna. The 
sparse data generated by SLA is first decomposed using 
IFKS A technique to obtain the signal subspace and later, the 
ESPRIT is applied to estimate the DOAs. Simulation results 
revealed that, for the ULA when compare to ESPRIT the 
proposed method performs better. The proposed method is 
evaluated for various SLA geometries and the algorithm 
succeeds in estimating the DOAs for all the cases. However, 
further investigation has to be made to handle the coherent 
signal sources. 
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Figure. 7. Histogram of the estimated DOAs fro case 2. 
Assumed DOAs of two signals [id 3 , 45 3 ] SNR = 10 dB 
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Figure. 8. Histogram of the estimated DOAs fro case 3. 
Assumed DOAs of two signals are [l[P. 45 11 ] SNR = 10 dB. 
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